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Gamma-ray emission from Wolf-Rayet binaries 
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Abstract. In the colliding wind region of early-type binaries, electrons can be accelerated up to relativistic ener- 
gies displaying power-law spectra, as demonstrated by the detection of non-thermal radio emission from several 
WR+OB systems. The particle acceleration region, located between the stars, is exposed to strong photon fields 
in such a way that inverse Compton cooling of the electrons could result in a substantial high-energy non-thermal 
flux. In particular cases, the ratio of the energy densities of magnetic to photon fields in the colliding wind region 
will determine whether a given source can produce or not significant gamma-ray emission. We present here a 
study of the binaries WR 140, WR 146, and WR 147 in the light of recent radio and gamma-ray observations. 
We show that with reasonable assumptions for the magnetic field strength WR 140 can produce the gamma-ray 
flux from the EGRET source 3EG J2022+4317. WR 146 and WR 147 are below the detection threshold, but 
forthcoming instruments like INTEGRAL and GLAST could detect non-thermal emission from them. 

Key words, stars: early-type - stars: binaries - stars: winds, outflows - radio continuum: stars - gamma-rays: 
observations - gamma-rays: theory 



1. Introduction 

Early-type stars present strong supersonic winds that are 
responsible for significant mass loss rates, which in Wolf- 
Rayet (WR) stars can reach values close to 10 -4 Mq yr _1 
(Abbott et al. |1986| , Leitherer et al. |1997| ). These winds 
interact with the interstellar medium sweeping up the am- 
bient material, creating cavities or bubbles (e.g. Benaglia 
& Cappa 1999| ), and forming strong shock fronts. 

At the terminal shocks, with typical velocities of thou- 
sands of km s _1 , locally-injected supra-thermal particles 
could be accelerated up to relativist ic ene rgies and power- 
law d istributions (e.g. Casse & Paul 1980] , Volk & Forman 
1982 ) . Energetic particles from the star lose too much en- 



at 



energies 
l 



E > 100 MeV is in the 



massive stars 

range i() 32-33 erg s _1 (Benaglia et al. |200l| ), too low 
to be detected at standard distances, by instruments 
like the Energetic Gamma-Ray Experiment Telescope 
(EGRET) of the recently terminated Compton Gamma- 
Ray Observatory mission or even by the forthcoming 



AGILE Italian satellite (Mereghetti et al. 2001) 



ergy in the expanding wind before reaching the shock to 
bo efficiently accelerated, but it has been suggested that 
partial re-acceleration during their travel, produced by 
multiple shocks from line-driven instabilities in the in- 
ner wind region, could compensate the adiabatic losses 
(White 1985). In any case, if a continuous injection of 
supra-thermal protons or ions can be sustained, the sub- 
sequent interaction of the relativistic hadrons with ambi- 
ent atoms will produce gamma-ray emission through the 
neutral pion chain: p + p — > 7r° + X , 7r° -^7 + 7. 

For typical densities in the interstellar medium, how- 
ever, the expected gamma-ray luminosity from isolated 



In the case of early-type binaries the situation is dif- 
ferent. In a typical WR+OB binary the collision of the 
winds from both stars produces a strong shock at some 
point between the stars, in a region exposed to strong 
UV stellar fields. Both electrons and protons can be ac- 
celerated in these colliding wind regions (Eichler & Usov 



1993 ). Electrons will cool through synchrotron and inverse 



Compton (IC) radiation. In fact, the detection of non- 
thermal radio emission in many early-type binaries cor- 
roborates the existence of a population of relativistic elec- 
trons in some of these systems (e.g. Dougherty & Williams 
|200C and references therein). In some cases, like Cygnus 
OB2 No. 5, the colliding wind region is spatially resolved 
with VLA observations and it appears as an extended, 
lightly elongated non-thermal feature in the radio images 
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(Contreras et al. 1997). In such cases, the geometry of the 
system can be established, radio flux components can be 
separated, and detailed calculations of the expected IC 
emission at high energies can be done. For the particu- 
lar case of Cyg OB2 No. 5, for instance, Benaglia et al. 



(2001) estimate that about a half of the gamma-ray flux 
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measured by EGRET from the Cyg OB2 association could 
come from the early-type binary. 

Using data from the third and final EGRET catalog 
of point-like gamma-ray sources (Hartman et al. 1999), a 
correlation analysis between unidentified sources on the 
one hand and WR stars on the other shows several co- 
incidences (Romero et al. 1999). The a priori chance 
probability is estimated in the range 1CP 2 — 10~ 3 , which 
is suggestive, but certainly not overwhelming as for the 
case of supernova remnants (Romero et al. 1999, Romero 
2001). Some positional associations, however, deserve fur- 
ther study in the light of the most recent observations. 
This is particularly true for WR 140, a very interesting 
system whose possible high-energy emission has been al- 
ready discussed in the pre-E GRET era by Pollock ( 1987 ) 
and also by Eichler & Usov ( 1993 ) in their seminal paper 
on non-thermal radiation from WR+OB binaries. The lo- 
cation of WR 140 is within the 98% confidence contour 
of the gamma-ray source 3EG J2022+4317. The existence 
of recent and detailed radio observations of this WR bi- 
nary provides new tools for a reliable estimate of the ex- 
pected gamma-ray flux, which can be now compared with 
the measured EGRET flux, not available at the time of 
Eichler & Usov's paper. 

Physical conditions in WR binary systems with non- 
thermal colliding wind regions can be probed through 
gamma-ray observations. When thermal and synchrotron 
components of the radio emission can be adequately sep- 
arated, the orbital parameters are well known, and the 
spectral types of the stars are determined, then the ratio 
of the synchrotron to the gamma-ray luminosity can be 
used to estimate the magnitude of the magnetic field in 
the shocked region. Even in case of a lack of clear gamma- 
ray detection, due to the instrument sensitivity, we can set 
bounds on the magnetic field strength, and hence make in- 
ferences on the field in the stars. The key point is that, 
granted the presence of a population of relativistic parti- 
cles and an UV-photon field, both of which we know are 
present in WR binaries with non-thermal radio emission, 
then the production of IC gamma-rays is unavoidable. The 
question, of course, is whether these gamma-rays come 
with a flux density high enough as to be detected by the 
current technology. 

In this paper we will study the case of three WR bi- 
naries: WR 140, WR 146, and WR 147. We will use the 
most recent results from radio observations in order to fix 
the value of model parameters, and then we will compute 
the expected gamma-ray luminosity with reasonable as- 
sumptions for the magnetic fields. We then will compare 
with gamma-ray observations in order to test our original 
assumptions. Finally, within the constrains imposed by 
the observations, we will make some predictions for future 
gamma-ray instruments like INTEGRAL and GLAST. 

We will not discuss cases of very close (i.e. short pe- 
riod) binary systems because they have typical separations 
of ~ 0.1 AU, which impose an upper limit to the size of the 
particle acceleration region at the colliding wind shocks 
that results, in turn, in a severe constraint for the high- 



est possible energy of the relativistic electrons. Even with 
strong magnetic fields of ~ 1 Gauss in the shocked re- 
gion, electrons could not go well within the GeV domain. 
However, soft gamma-rays of a few MeV are still possible, 
and some of these systems should deserve further high- 
energy studies as potential targets for the INTEGRAL 
satellite, which is optimized for such an energy range. 

The structure of the paper is as follows. In the next 
section we will present the model for the gamma-ray emis- 
sion from the colliding wind region. In Section ^ we will 
discuss the main characteristics of the WR binaries in our 
sample. Section |^ presents the existing gamma-ray infor- 
mation about the regions where the stars are located. In 
Section || we give the main results whereas in Section |^ 
we discuss some of their implications. We close in Section 
with the conclusions. 

2. Gamma-ray emission from massive binaries 

Gamma-ray production in early-type binaries with collid- 
ing winds has b een d iscussed by Eichl er fc U sov ( 1993 ), 
White fc Che n ( |1995| ), Benaglia et al. fl200l| ) and Miicke 
& Pohl ( 2002 ). We will present here the main features of 
the basic model. 

In an early-type binary system the winds from the pri- 
mary (e.g. a WR) and the secondary (e.g. an OB) stars 
flow nearly radially and collide at a point located at a 
distance ri from the i-star, given by : 
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In these expressions the subscript "1" stands for the pri- 
mary star, and "2" for the secondary, D is the binary 
separation, and the parameter r/ is defined in terms of the 
wind terminal velocities Uqo and the stellar mass loss rates 
M: 



M?v n 



V 



Mi v .. 



(2) 



1^00,1 



At the colliding wind region, first-order diffusive shock 
acceleration (e.g. Drury 1983 ) results in the production 
of a power-law spectrum AT (7) oc 7 _p , with p ~ 2 for 
electrons of energy E = jm c c 2 . Because of the strong 
photon fields existing in the colliding wind zone, IC losses 
are expected to dominate the radiative cooling at high 
energies. The high energy cutoff of the electron spectrum 
can be then obtained from the condition that the IC loss 
rate for the particles does not exceed the acceleration rate 
(Eichler & Usov |l993| ): 
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Here B is the local value of the magnetic field and 
Lob is the bolometric luminosity of the OB star, which 
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is closer to the shock. Typical values for the maximum 
Lorentz factor are, according to Eq. (||), of a few times 
10 4 . 

Inverse Compton losses will also produce a modifica- 
tion in the spectrum of the relativistic particles. A break 
is expected at the energy at which the cooling and escape 
times are equal (e.g. Longair 1997, p. 281). This will occur 
at a Lorentz factor 7b given by: 



7b 



3m e c 
Aa T Ut c 



(4) 



where U is the average energy density of the photon field, 
t csc is average time spent by the particles in the field re- 
gion (typically t csc ~ s /c, with s the linear size of the col- 
liding wind region), m c is the electron mass, and <tt the 
Thomson cross section. The spectrum will steepen from 
an index p to p + 1 for energies higher than 7^,. 

The IC photons produced in a stellar photon field 
with seed photons of frequency z/» = 5.8 10 10 T e g, where 
T e ff is the stellar effective temperature, are characterized 
by v lG = 4/37 2 ^*. The electrons will also lose energy 
through synchrotron emission in the magnetic field, pro- 
ducing a synchrotron flux with a spectrum S{y) oc v^ a 1 
with a — (p — l)/2. The frequency of the synchrotron 
photons will be v syn = 4.2 B -f 2 MHz. Depending on the 
strength of the magnetic field, the tail of the synchrotron 
emission can reach even optical frequencies in some cases. 
Of course, in such bands it is not observable because of 
the overwhelming thermal contribution, but at radio wave- 
lengths it can be detected and measured with interfero- 
metric instruments. 

The rate of IC interactions in the colliding wind region 
per final photon energy is: 
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where eo is the typical energy of the seed photons. For the 
case of Thomson scattering and adopting the "head-on" 
approximation (i.e. the photons- an isotropic field in the 
lab frame- are treated as coming from the opposite direc- 
tion t o the electron's velocity in the electron frame, see 
Jones 1968 for details), the function f{x) can be written 
as: 



f(x) = -(l-x)P(l/4 7 2 ,l,x), x 



4e 7 2 



(6) 



The specific luminosity is then obtained by integrating 
the scattering rate in Eq. (j^) over the particle energy dis- 
tribution, and multiplying by the observed photon energy 
e and the number density n p h = U/e$. This leads to the 
expression: 

dL kV a^m^f? U e 
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Here, k is the constant involved in the energy distri- 
bution of the electrons, 72 the highest energy of the elec- 
trons considered in the calculation, and V the interacting 
volume. According to the energy range we are interested 
in, 72 can be equal to (or smaller than) 7 max . When p 
is steeper than 1 (as in our case) and 7 m j n -C 72, this 
expression reduces to: 



L f = 



dL kVa T m e c 3 U2P- 
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A direct calculation of the IC luminosity of the col- 
liding wind region of a particular binary system requires 
a knowledge of the parameter k in the electron energy 
distribution. This parameter is unknown, but we can use 
information on the synchrotron luminosity of these same 
electrons in order to circumvent this problem. In partic- 
ular, the ratio of synchrotron to IC luminosities for the 
colliding wind region can be estimated as (White & Chen 



1995, Benaglia et al. 2001): 
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where B is expressed in G, is the distance to the i-star 
in AU, and Li is the star luminosity in Lq units. Since 
L S ym T"i, and Li can all be, at least in principle, determined 
from observations, the gamma-ray visibility of the massive 
binary will crucially depend on the value of the magnetic 
field in the particle acceleration region. 

It is usually assumed that the external magnetic field 
of the star in the absence of stellar wind is dipolar, and 
that with wind it obeys the standard B oc radial de- 
pendence for large r given by Eichler & Usov (1993): 



B ~ £L 



valid for r > R* — — . 

rxr V to t 



(10) 



where i?* is the star radius, £?* the surface magnetic field, 
K-ot the surface rotational velocity, and rx ~ R* the 
Alfven radius. The value of the surface magnetic field is 
not well known, but it is thought that in WR stars could 



reach 10 4 G (e.g. Ignace et al. 1998). Then, in the colliding 
wind region of typical WR+OB binaries, the field could 
be in the range 10~ 2 — 10 G. 

At the colliding wind zone an eq uipar tition magnetic 
field can be derived following Miley ( 1980 ): 
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Here 6 X , 9 y represent the source sizes in mas; x is t ne 
energy ratio between heavy particles and electrons; / is 
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the filling factor of the emitting region; s is the path length 
through the source in the line of sight, in AU; (j> is the angle 
between the uniform magnetic field and the line of sight; 
F the flux density in mJy of the region at frequency v§\ 
and vi and v 2 are the upper and lower cut off frequencies 
presumed for the radio spectrum, these last three in GHz. 
For minimum energy conditions, X = / = sin</> = 1. 

IC scattering of stellar photons is not the only mech- 
anism capable of producing gamma-rays in the collid- 
ing wind region of an early-type binary. Relativistic 
bremsstrahlung in the ions of the winds will also result into 
gamma-ray emission. On the other hand, the same first- 
order acceleration mechanism that forms the electronic 
relativistic population should operate upon protons and 
ions. The decay of neutral pions generated in the hadronic 
interactions between the relativistic protons and the nu- 
clei in the wind yields gamma-rays, which at high energies 
present the same energy spectrum than the parent pro- 
ton population. At 67.5 MeV the spectrum should present 
the typical pion bump. All these complementary contribu- 
tions, however, are rather minor in comparison with the IC 
emission because of the ambient densities involved in col- 
liding wind binaries. In Section ||, when we will present the 
results of our calculations of the IC gamma-ray luminos- 
ity for several systems, we will also provide the expected 
bremsstrahlung and hadronic gamma-ray luminosities for 
completeness, but the reader is referred to Benaglia et al. 
(2001) and references therein for the formulae. 



3. Selected binaries 

The recently updated catalogue of WR stars (van der 
Hucht 2001) lists about 230 objects, all but one at dis- 



tances up to 20 kpc. Two main radio surveys with high- 
resolution have bee n carr ied out. The first one, using the 
VLA (Abbott et al. 1986 ) compiled data from ~ 40 stars 
up to 3 kpc, north of declination —47°. The southern one 



wa s don e using ATCA (Leitherer et al. 1997 , Chapman et 
al. |1999| ) and encompasses ~ 40 WR stars with 5 < 0° and 
distances up to 3 kpc. Special obje cts w ere re-observed in 
recent year s (e.g. Contreras et al. 1997 , Dougherty et al. 



1996 , 2000 , etc). Some objects show spectral indices a ap- 
proaching the canonical thermal value of +0.6, whereas 
the rest present negative and/or variable values. 

For the present investigation, we looked first for candi- 
dates with detected non-thermal emission, indicating the 
presence of relativistic electrons. Taking into account the 
works mentioned above, a list can be made containing the 
objects that have shown spectral indices other than ther- 
mal. Such a list is presented in Table ^, where we have 
separated the stars in three groups or types, according to 
the spectral index: 'Comb' sources are those that present 
an index near or variable; 'NT' ones have a clearly nega- 
tive value of a. Two systems - WR 146 and WR 147- have 
been 'resolved' using VLA and MERLIN in two sources 
each: one thermal, the other not. The spectral classifica- 
tion, distances, and periods in the table are taken mostly 



the spectral indices. We also provide information on the 
different wavelengths at which the stars were detected, 
whenever available. 

After eliminating objects with no evidence of colliding 
winds or close binaries, we are left with WR 146 and WR 
147 as the best candidates. The colliding wind region has 
been resolved in continuum radio observations for these 
systems. They have orbital periods of hundreds of years. 
The geometry of the colliding wind region is relatively 
clear for these binaries. 

A third classical candidate was added to our sample: 
WR 140. Although the colliding wind region is not well- 
resolved here, this star has been monitored over its com- 
plete period, it has been pointed as the probable counter- 
pa rt of an EGRET gamma-ray point source (Romero et 
al. 199E ), and it wa s also the original example discussed 
by Eichler & Usov ( 1993 ) to illustrate their non-thermal 
model. 



3.1. WR 140 

WR 140 (HD 193793, V1687 Cyg) is a spectroscopic bi- 
nary system, formed by a WC7pd and an 04-5 compan- 
ion (van der Hucht 2001 ). It was the first star for which 
non-thermal emis sion has been detected (Florkowski & 
Gottestman 1977 ). The condition o f colli ding wind binary 
was confirmed by Williams et al. ( |1987| ), deducing a pe- 
riod of 7.9 yr and an eccentricity e = 0.84. It has shown 
variability at radio, IR, optical, UV, and X-rays (Setia 
Gunawan et al. 2001b and references therein). Periodic 
dust formation around periastron passage would be re- 
sponsible for the IR excess observed (Williams et al. |1987D . 



Both White & Becker (1995) and Williams et al. (|1990| 



1994) monitored the system at radio frequencies. The first 



from van der Hucht ( 2001 ) . Special references are given for 



authors used the VLA at 2, 6 and 20 cm, from 1985 to 
1993. Williams et al. observed WR 140 from 1986 to 1994 
with WSRT. The observations allo w to e stimate a spectral 
index of —0.6 (van der Hucht et al. 2001 ). They adopted a 
model in which the WR stellar wind is strongly enhanced 
at the equatorial plane, where most of the mass loss is con- 
fined. The emission suffers from varying circumstellar ex- 
tinction in the line of sight to the non-thermal radio source 
that is identified with the colliding wind region. The orbit 
s hould be inclined. More recently, Setia Gunawan et al. 
( 2001b| ) presented the results of 16 years of UV monitor- 
ing of WR 140, and provided explanations for the observed 
spectral variability. 

The stellar parameters of WR 140 relevant to this work 
are listed in Table ^| along with the corresponding refer- 
ences. Throughout the period, the separation between the 
stellar components of the system ranges from 2.4 to 27 
AU. We have assumed a separation D = 10 AU in or- 
der to compare with gamma-ray observations (see Section 
^), and computed a value r\ = 0.06, that was used in our 
calculations (Section |)|). The WR stellar luminosity was 
taken as an average of all model luminosities given by 
Koesterke & Hamann ( |1995 ) for WC7 stars. The luminos- 
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Table 1. Wolf-Rayet stars with non-thermal -or combined- emission. 



WR 


Sp. Class. (a) 


d 

(kpc) 


Period (a) 
(days) 




Type 


Observations 
(cm) 


1 1 
1 1 


Wjr^o i r\n trTTT \t 
WW + Lm.OIII-V 




7Q r 9 

lo.oo 


[-0.5J* 


Comb 


9 a i 9 on 
o,0,lo,zU 


14 


WC7 + ? 


2.00^ ' 


2.42 


1 n(b) 
— 1.0' ; 


NT 


3,6 


on 

39 


WUI + (Jti! 


0.t>3 v ' 




n n(b) 


Comb 


9 t; i 9 on 
3,6,13,20 


48 


WC6(+O9.5/B01ab) 


2.27 v ' 


18.34 


—0.4, [ — U.8J ^ ' 


IN 1 


3,6,13 


nn 
90 








n nfb) 
0.0 1 ' 


Comb 


3,b,13 


105 


WlN9n 


1.58^ ' 




n o f b) 

—0.3' ' 


Oomb 


3,6 


112 


WC9d + UB : 


/lie (a) 




i ocrfb) 

— 1.35' ' 


IN 1 


3,13 


125 


WC7ed + 09111 


3.06 ( - a - ) 


>6600 


—0.5 — > 0.7'' ' 


Comb 


2,6,20 


137 


WC7pd + 09 


2.38 (a) 


4765 


0.0 (c) 


Comb 




140 


WC7pd + 04-5 


1.10 (a) ,1.3 (f) 


2900 


-0.5 (g) -0.6 (h) 


Comb 


2,6,20 


146 


WC6 + 08 


0.72 (a) ,1.25 (h) 


ELPB (i) 


-0.62 (j) 


RESOLVED 


1.3,3,6,20 


147 


WN8(h) + B0.5V 


0.65 (a) 


ELPB (i) 


-0.43 (k) ,-l.l (1) 


RESOLVED 


1.3,2,3.6,20 



(j) Dougherty et al. 2000|; (k) Setia Gunawan ct al. 2001a; (1) Skinner et al. 1999 



1995c (h) van der Hucht et al. 2001 ; (i) stands for Extr emely Long Period Binary; 



(*) Brackets mea n non -thermal index deri ved fr om model; (a) van der Hucht 2001 ; (b) Chapman et al . 1999; 
(c) Abbott et al. [l986|; (d ) Williams et al. |l992| ; (e ) Dou gherty & Williams |200C| ; (f) Smith et al. [l990| ; 
(g) White & Becker 



ity class of the secondary is suggested by Setia Gunawan 
et al. ( |2001b| ) as 04-5 (V). The effective temperature and 
radius of the primary were assumed the same of WR 125 



(Koesterke & Hamann 1995), because of their similarities. 
The effective temperature, stellar luminosity, and radius 
for the se condary were taken from the tables of Vacca et 
al. ( 1996| ), for a spectral type 04.5 V. A mean molecular 
weight (/i) of 5 was adopted for the primary, and 1.5 for 
the secondary. The mass loss rate of the sec ondar y was 
computed following the models of Vink et al. (2000). 



Van der Hucht et al. (2001) give values for mass loss 
rates and wind terminal velocities of both components, 
the spectral non-thermal index, distance, and non-thermal 
flux density at 6 cm. In order to obtain an estimate of the 
extension of the colliding wind region, we have computed 
the radii of the photospheres at 5 GHz of both compo- 
nents, using the expression from Wright & Barlow ( |1975| ). 
After verifying they overlap, we have considered the size 
of the non-thermal source equal to the OB radio photo- 
sphere, about 4 AU. 

The equipartition magnetic field can be calculated 
with Eq (pT|) using v x = 1.4 GHz and v 2 = 250 GHz 



source with 6 X = 6 y = 3 mas, and s = 4 AU. Inaccuracies 
in the colliding wind region size or in the filling factor in- 
troduce the main errors in the determination of the mag- 
netic field. For example, if / = 0.1, or s = 2 AU, the 
magnetic field doubles its value. 

3.2. WR 146 

This system contains the brightest WR star detected in 
radio continuum so far. The secondary is an OB star in 
an orbit with a period estimated in ~ 300 yr (Dougherty 



et al. 1996) 



Afte r obs erving with the EVN and VLA, Felli & 
Massi (1991) detected emission with a spectral index 
a = — 1. Dougherty et al. (1996) resolved the system us- 



ing MERLIN in two sources: a bright non-thermal north- 
ern component of 38 ± 1 mas in diameter, and a weaker 
southern component. Niemela et al. (1998) reported the 



resolution of the system after HST-WFPC2 observations. 
Their result for the binary separation was confirmed by 22 
GHz observations taken by Dougherty et al. (2000), who 



(see Wendker 1995), and assuming a spherical non-thermal 



detected the stellar wind of the companion star with the 
VLA. Their radio maps allowed to derive a value for r 2 
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Table 2. Stellar parameters adopted for the WR 140 system 



Variable 


WR 




OB 


unit 


Separation from colliding wind region 


8.5W 




1.5<* 




AU 


Spectral Class. 


WC7pd (a) 




04-5 (V) (b) 




log (L/L ) 


5.35 (c) 




5.80 (d) 




T c ff 


48000 (c) 




47400 (d) 


K 


R, 


6.0 (c) 




12.0 (d) 


1' 


/< 


5(0) 




1.5 (c 








2860 (c) 




3100 (c 




km s _1 


M 


5.7 10 _5(e) 




3.0 10" 6(c 




M yr- 1 


Distance d 




1320 te) 






pc 


V 




0.06 








Non-thermal spectral index 




-0.60 (e) 








SsGHz 




22.5 (c) 






mJy 


Size of non-thermal source 




6.0 10 13(c) 






cm 


Photon energy density U 




38.5 






erg cm -3 


Bcwr (equipartition) 




200 (c) 






mG 


(*) for binary separation of 10 AU; (a) 


van der Hucht 2001; (b) Set 


ia Gunawan et al. 2001b; 


(c) adopted here, see text; (d) Vacca et al. 1996 ; (e) 


van der Hucht et al. 2001 







of ~ 50 AU, at d = 1.25 kpc, and 77 = 0.1. They took 
optical spectra of the system and interpreted the results 
as strongly suggesting that the 08 star is an Of, possibly 
a supergiant. 

Setia Gunawan et al. ( 2000| ) presented the results of 
1.4 and 5 GHz observations with WSRT, from 1989 to 
1999. They found three different kinds of flux variabil- 
ity behaviour on this time-span: a linear increasing trend 
observed during 10 years, a possible 3.38-yr periodic sig- 
nal, and rapid fluctuations on time scales of weeks; they 
also derived a time-averaged non-thermal spectral index 
a = -0.62. 

The values adopted for our calculations are shown 
m Table § Since the distance estimates to the WR 146 
binary system differ from 0.75 kpc to 1.7 kpc (Setia 
Gunawan et al. |2000 and references therein), we have car- 
ried out the calculations ad optin g the mean value d = 1.25 
kpc (van der Hucht et al. 2001 ). The corresponding sep- 
aration of the two stars is 210 AU (Setia Gunawan et al. 
2000| ) . The stellar luminosity of the secondary is an aver- 
age of the values given by Setia Gunawan et al. ( [2000 ). The 
effective temperature and radius of th e seco ndary were 
taken from the tables of Vacca et al. ( |1996[ ). The mass 
loss rate of the secondary was computed from Eq. (§). 
The colliding wind region size was derived assuming an 
extension of the non-thermal source of 46 AU at 1.25 kpc 



(Setia Gunawan et al. |2000| ). We estimated the equiparti- 
tion magnetic field between 327 MHz and 22 GHz (Taylor 
et al. [1996D, with 6 X = 9 y = 38 mas, and s = 47 AU. 



3.3. WR 147 

Non-thermal emission form this system was reported after 



VLA observations by Abbott et al. (1986). Higher angu- 
lar resolution observations taken with the same instru- 



servations by Moran et al. (1989) confirmed the presence 
of two sources: a southern thermal one (WR 147S) su- 
perposed with the WR star, and a northern non-thermal 
component (WR 147N). Williams et al. ( |1997[ ) observed 
the system again with MERLIN and found that the two 
sources are separated by 575 ± 15 mas. They reported 
also a faint IR source near WR 147N, but slightly far- 
ther to the WR star than the radio source WR 147N, and 
derived a spectral type B0. 5 V for the IR star (con firme d 
later through optical observations by Niemela et al. 1998 ). 
Because of the presence of non-thermal emission between 
the two stars and its location much closer to that with the 



weaker wind, Williams et al. ( 1997 ) proposed the system 
as a colliding wind binary. 

Using VLA observations at 3.6 cm, Contreras & 
Rodriguez ( 1999| ) found that the non-thermal wind- 
interaction zone remained constant in flux density during 
1995-1996. The thermal emission from the WR star, on 
the contrary, increased ~ 25%, probably reflecting the in- 
homogeneous nature of the wind. Both the observed radio 
morphology and the theoretical modeling by Contreras & 
Rodriguez ( 1999] ) clearly supports a colliding winds sce- 
nario for WR 147. 



ment (Churchwell et al. 1992) and MERLIN 5-GHz ob 



Setia Gunawan et al. ( 2001a ) presented the results of 
a monitoring campaign of the source from 1988 to 1997 at 
1.4 and 5 GHz with WSRT. Once subtracted the south- 
ern thermal contribution, the spectral energy distribution 
could be fitted by a synchrotron emission model which 
includes free-free absorption. Flux density variations on 
different time scales can be explained by considering inho- 
mogeneities in the wind, plasma outflow, etc. After fitting 
the thermal source data, they computed a non-thermal 
flux density Snt = Stotai — >Sfitj obtaining the spectral in- 
dices ^thermal = +0.6 and a syn = —0.43 for the thermal 
and synchrotron components, respectively. The statistical 
significance of this result, however, is very low. Skinner 
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Table 3. Stellar parameters adopted for the WR 146 system 



Variable 


WR 




OB 


unit 


Separation from colliding wind region 


160 W 




50 (,) 


AU 


Spectral Class. 


WC6 (a) 




08 If (b) 




log (L/L Q ) 


5.10 (c) 




5.00 (d) 




T e g 


49000 (o) 




35700 (f) 


K 


H, 


5.0 (c) 




23.1< f » 


R© 


/' 


5.29 (s) 




1.34 (h) 






2700 (g) 




1300 (i) 


km s" 1 


M 


2.6 10 _5(g) 




5.4 10" 6(d) 


M yr" 1 


Distance d 




1250 w 




pc 


V 




0.10 






Non-thermal spectral index 




-0.62 (c) 






SsGHz 




28.5 (i) 




mJy 


Size of non-thermal source 




7.0 10 14 (d) 




cm 


Photon energy density U 




6.1 10" 3 




erg cm -3 


B cwr (equipartition) 




25.0 (d) 




mG 



(*) deduced for r 2 = 40 ± 9 mas (Dougherty et al. |200C|); (a) Smith et al. |199(| : 



2000 



(b) Dougherty et al. 
(e) Koesterke & Hamann 



(c) S etia Gunawan e t al 
(f ) Vacca et al 



2000 



1996 



(d) adopted here, see text; 



fg) W illis et al. 



(h) Lamers & Leitherer 1993; (i) van der Hucht et al. 2001 



1997 



et al. (1999) have concluded, using nearly simultaneous 
observations at five frequencies, that there is a significant 
steepening in the non-thermal spectrum above 10 GHz. 
These authors have even suggested that a monoenergetic 
relativistic electron spectrum injected at the source might 
be a plausible explanation of the observed radio flux dis- 
tribution. The quality of the radio data, unfortunately, 
does not allow to draw any firm conclusion. 

The parameters for the system WR 147 adopted in our 
calculations are listed in Table ^, with the respective ref- 
erences. A canonical value of a = —0.5 was adopted, cor- 
responding to the expected p = 2 spectrum in the electron 
distribution at relatively low energies. Compton losses will 
introduce a steepening at higher energies. In Section || 
we will revisit this assumption and we will briefly discuss 
the implications of the monoenergetic electron distribu- 



tion suggested by Skinner et al. (1999) for the high-energy 



emission. The equipartition field at the colliding wind re- 
gion was computed for 6 X = y = 267 mas, and s = 175 
kpc, b etween frequencies 350 MHz (S ctia Gunawan et al. 
2001a| ) and 250 GHz ( Wendker [l995 ) . 



4. Gamma-ray observations 

The gamma-ray observations used here were provided by 
the EGRET telescope of the Compton satellite. The re- 
duced data were published in the third EGRET catalog 



(Hartman et al. 1999). The catalog lists the best estimated 
position of 271 point gamma-ray sources, along with the 
location error boxes (taken as the 95% confidence con- 
tours), integrated gamma-ray fluxes in the energy range 
100 MeV - 20 GeV, spectral indices, and other informa- 
tion. 



Only WR 140 is positionally coincident with a point- 
like source in the catalog, 3EG J2022+4317. The pho- 
ton flux summed over satellite observing cycles 1, 2, 3, 
and 4 is 24.7 ± 5.2 10 -8 ph cm -2 s -1 . The photon spec- 
tral index F(E) cx is T = 2.31 ± 0.19. The loca- 
tion error box is large, with a radius of ~ 0.7 degrees, 
which renders a clear identification almost impossible in 
practice. However, WR 140 is the only highly energetic 
source known in this region (Romero et al. |l99S|) . Figure 
shows the relative position of WR 140 with respect to the 
gamma-ray source. 

The EGRET viewing periods for this source cover since 
1991 to 1994. Following the orbital positions and times 



deduced by Setia Gunawan et al. (2001b), the stars were 
separated between ~ 4 to ~ 16 AU during the time-span 
1991.5 - 1994.5. We took a mean distance D = 10 AU 
for the binary separation during the gamma-ray observa- 
tions in order to perform the calculations of the expected 
average flux. 

The variability analysis of the gamma-ray data for 
3EG J2022+ 4317 g ives a variability index I = 1.3 in 
Torres et al. ( [2001a| ) scale, which is normalized to the av- 
erage pulsar variability. A source will be variable if / > 5, 
at 8cr level. Non-varia ble sou rces have in dices typically 
I < 1.7 (Torres et al. |2001b| ). Tompkins ( |1999| ) analysis 
yields for this source values of his r index of < r >= 0.13, 
Tmin = 0.00, and T max = 0.50. Here, a variable source will 
have T m ; n > 0.6 (see Torres et al. 2001b for a comparison 



of both indices). Consequently, the source has not been 
formally variable during the observing time interval. But 
a note of caution is necessary here. As it can be seen in 
Fig. ^, upper panel, all but 3 observing periods for this 
source yielded only upper limits. In such a case the for- 
mal variability index can be misleading because under the 
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Table 4. Stellar parameters adopted for the WR 147 system 



Variable 


WR 




OB 


unit 


Separation from colliding wind region 


374(a) 




43(b) 


AU 


Spectral Class. 


WN8(h) (c) 




B0.5 V (a) 




log (L/L ) 


5.52 (d) 




4.70 (c) 




T e g 


26000 (d) 




28500 (f) 


K 


R>. 


20.6 (e) 




8.0 (g) 


I' 


/' 


3.18 (c) 




1.5< b » 






950 (e) 




800 (b) 


km s _1 


M 


2.4 10" 5(e) 




4.0 10- 7 < b) 


M yr- 1 


Distance d 




650 (c) 




pc 


V 




0.014 






Non-thermal spectral index 




-0.5« 






SsGHz 




12.5< h > 




mjy 


Size of non-thermal source 




2.6 10 15(b) 




cm 


Photon energy density U 




4.0 10" 3 




erg cm -3 


-Bcwr (equipartition) 




5.0« 




mG 



(a) Williams et al. 1997; (b) Setia Gunawan et al. [2001a , see text for a discussion; (c) Smith 



et al. 



!!)<)(; 



al. 1996; (h) van der Hucht et al. 2001; (i) adopted here, see text 



(d) Crowther et al. 1995 ^(e) Morris et al. 2000; (f) Crowther 1997; (g) Vacca et 



I 

a 
-I 



O 




82.0 81.0 80.0 

Galactic Longitude 

Fig. 1. EGRET probability contours for 3EG 
J2022+4317. Contour labels are 50% 68%, 95% and 
98%. The position of WR 140 is marked 



upper limits real variability might be hidden. In the lower 
panels of Fig. || we show the radio evolution of the WR 140 
for the same epoch of the gamma-ray observations. The 



data is taken from White & Becker (1995), directly from 
1991 to 1993.6 and by extrapolation from 1993.6 to 1994.5. 
We can see that the radio flux was significantly decreas- 
ing from 1992 to 1993. If the value of the gamma-ray flux 
of the source was close to the upper limit established by 
EGRET in mid- 1992, then the high-energy emission might 
also have decreased significantly that year. The poor res- 



olution of the EGRET lightcurve prevents any conclusion 
in this respect. It is not possible, in our opinion, to reject 
the association of WR 140 with 3EG J2022+4317 on the 
sole basis on the low variability indices. 

Although the authors of the third EGRET catalog give 
the approximate upper limits for gamma-ray sources at al- 
most any point in the sky, at the positions of WR 146 and 
WR 147 no threshold could be computed. The reason is 
that the positions of the stars are close to the fairly strong 
source 3EG J2033+4118. This makes it impossible for the 
gamma reduction programmes to get a statistically mean- 
ingful upper limit for weaker sources in the surroundings. 
Because of this contamination problem, we consider for 
the positions of WR 146 and WR 147, an EGRET detec- 
tion limit of 70% of the flux of the gamma-ray source 3EG 
2033+4118, i.e. 50 10~ 8 ph cm -2 s~\ which seems to be 
a reasonable assumption (Bob Hartman, private commu- 
nication). 

5. Results 

The following results for the different WR binaries in our 
sample are obtained with the parameters listed in Tables 
§ - 0. A summary of the results is presented in Table |[ 
where we list, for each system, the expected gamma-ray 
luminosity due to inverse Compton scattering, the rela- 
tivistic bremsstrahlung, and the pion-decay contributions 
to the gamma-ray flux generated at the colliding wind re- 
gion. In the following subsections we describe the main 
results of our calculations for each star. 



5.1. WR 140 

For this system we have -B cl 
photon energy density of U 



~ 0.2 G and a very high 
38.5 erg s _1 . The inverse 
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Fig. 2. Light curves at a) gamma-rays, b) 2cm, and c) 6cm. 



Compton losses impose an upper energy limit for the lo- 
cally accelerated electrons of 7 max = 4.2 10 4 and a break 
in the spectrum at 7b ~ 400. The break will appear in 
the synchrotron spectrum at syn w 135.2 GHz and at 
gamma-ray energies at = 2 Ah MeV, below the EGRET 
energy range. The Lorentz factors of electrons radiating 
IC gamma-rays in the ambient photon fields and con- 
tributing to 100 MeV-20 GeV are between 71 = 2.5 10 3 
and 72 = 3.6 10 4 . Electrons having Lorentz factors in 
this range are also capable of emitting synchrotron radi- 
ation between v x = 5.5 10 3 GHz and v 2 = 1.1 10 6 GHz. 
The total synchrotron luminosity from these electrons is 
- 2.5 10 30 erg s -1 . 



The total gamma-ray luminosity due to IC scatter- 
ing at the colliding wind region between 100 MeV and 20 
GeV results ~ 2.1 10 34 erg s _1 . At this energies, after 
the break in the gamma-ray spectrum, the photon spec- 
tral index should be T ~ 2.1 (original injection spectrum 
with p = 2.2 according to radio observations below the 
break) . Such an index is in reasonable agreement with the 
index observed by EGRET in the source 3EG J2022+4317 
(T = 2.31 ± 0.19), especially if we take into account that 
there are also uncertainties in the determination of the 
radio spectral index. 

Both the luminosities due to pion decays from 
p — p interactions (< 10 22 erg s^ 1 ) and relativistic 
bremsstrahlung in the winds (< 10 32 erg s -1 ) can be dis- 
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36,5 



36,0 



35.5 



34,0 



33,5 




Fig. 3. Gamma-ray luminosity (in erg s -1 ) between 100 
MeV and 20 GeV versus the magnetic field for WR 140. 
The horizontal line indicates the expected luminosity of 
3EG J2022+4317, if it is located at the distance of WR 
140. 



regarded in comparison to the IC luminosity (see Benaglia 
et al. 2001 for details of calculation). 

If the unidentified gamma-ray source 3EG J2022+4317 
is at the same distance of WR 140, the measured EGRET 



flux implies a luminosity of ~ 3.2 10 erg s , of the same 
order of magnitude than the computed luminosity due to 
IC scattering at the colliding wind region. Errors in the ob- 
served gamma-ray emission are at the level of ~ 25%. An 
additional factor that can affect the computations is the 
luminosity of the secondary star. For the WR 140 system it 
is not derived directly from observations, but interpolated 



from the tables of Vacca et al. (1996) and can be slightly 



overestimated. In any case, it seems likely that significant 
part of the gamma-ray flux of 3EG J2022+4317 might be 
contributed by the colliding wind region of WR 140. 

In Fig. |^ we present a plot of the gamma-ray luminos- 
ity of the colliding wind region of WR 140 versus the as- 
sumed magnetic field. The luminosity of 3EG J2022+4317 
at the same distance is indicated with a horizontal line. 
For high values of the magnetic field the synchrotron losses 
dominate over IC losses, and the gamma-ray emission is 
quenched. We emphasize, however, that in order to explain 
the gamma-ray source 3EG J2022+4317 through WR 140, 
no extreme hypothesis is required, and the Occam's razor 
principle is fulfilled since no new and otherwise yet un- 
detected object is postulated. In this sense, we suggest 
that WR 140 should be considered as the best currently 
available explanation for the origin of 3EG J2022+4317. 
The GLAST telescope, with its improved source location 
accuracy, will be able to test this proposition. 

5.2. WR 146 

For this system the equipartition magnetic field is signifi- 
cantly lower than in the case of WR 140: B CWI ~ 25 mG. 



35,0 



34,0 




B cwr (G) 

Fig. 4. Gamma-ray luminosity (in erg s _1 ) between 100 
MeV and 20 GeV versus the magnetic field for WR 146. 



This lower value is mainly due to the different size and 
location of the non-thermal regions. With the values of 
the parameters listed in Table [| we obtain the following 
results. 

The maximum energy of the accelerated particles in 
the colliding wind region is 7 max = 6 10 5 . These are very 
energetic particles indeed. At such energies Klein-Nishina 
effects are important. The highest energy for IC photons 
produced by these electrons is ~ 7 max m c c 2 /2 ~ 0.2 TeV. 

The Lorentz factors of the electrons contributing to 
EGRET's energy range are between 71 = 3 10 3 and 72 = 
4.2 10 4 . The local energy density U = 6.1 10~ 3 erg cm~ 3 
implies a break in the spectrum at very high energies, 
within the optical region (4.9 10 6 GHz). Electrons having 
Lorentz factors in the above range radiate synchrotron 
photons between v\ = 9.2 10 2 GHz and v 2 — 1.8 10 5 
GHz. The total synchrotron luminosity in this frequency 
range is 3.3 10 31 erg s — 1 . 

The gamma-ray luminosity due to IC scattering at 
the colliding wind region, contributed by seed photons 
from both stars is ~ 2.8 1 33 erg s _1 . In comparison, the 
gamma-ray luminosity due to pion decay is ~ 2.5 10 21 erg 



S _1 , and to relativistic brcmsstrahlung, ~ 4.6 1 29 erg s _1 , 
resulting both negligible. The IC spectral break appears 
at high energies (538 GeV), well above EGRET's range, 
so the spectrum in the MeV-GeV band should be much 
harder than in the case of WR 140, with values T ~ 1.6. 

Figure |] presents a plot of the gamma-ray luminos- 
ity between 100 MeV and 20 GeV versus the magnetic 
field. Very weak fields are ruled out by the EGRET non- 
detection. 

As we mentioned, the EGRET threshold towards the 
position of WR 146 can be taken as 50 10 -8 ph cm~ 2 s — 1 . 
This means that a gamma-ray source at the distance of 
the binary system would be detectable by EGRET if its 
luminosity is greater than 8.0 1 34 erg s _1 , which is more 
than an order of magnitude above the expected luminos- 
ity coming from the colliding wind region. This fact can 
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Table 6. Stellar magnetic fields 



0,000 



0,004 0,006 
Bcwr (G) 



Fig. 5. Gamma-ray luminosity (in erg s 1 ) between 100 
MeV and 20 GeV versus the magnetic field for WR 147. 



explain why no gamma-ray source was detected towards 
WR 146. However, as we will see, the expected flux is high 
enough as to be detected by GLAST, if the contamination 
problem from nearby sources can be solved. 

5.3. WR 147 

For WR 147 we have B cm ~ 5 mG. The maximum Lorentz 
factor for the electrons imposed by the IC losses is 7 max = 

1.0 10 5 . The high-energy cutoff frequency is then v C utoS = 

2.1 10 25 Hz or - 90 GeV. 

The Lorentz factors of the electrons whose IC emission 
falls in EGRET's energy range are between 71 = 3.3 10 3 
and 72 = 4.7 10 4 . Their synchrotron radiation is between 
v x = 230 GHz and v 2 = 4.6 10 4 GHz. The total syn- 
chrotron luminosity with such frequency interval results 
~ 5.7 1 30 erg s _1 . The break produced by the local pho- 
ton field of density U = 4.0 10 -3 erg cm -3 in the electron 
spectrum occurs at a Lorentz factor of 7b ~ 8.9 10 4 , and 
results in breaks in the synchrotron (at ~ 1.7 10 5 GHz) 
and IC spectrum (at ~ 72 GeV). 

The gamma-ray luminosity due to IC scattering at 
the colliding wind region, contributed by both stars, is 
~ 8.0 1 33 erg s _1 . In Fig. |^ we show how this luminos- 
ity would change with the assumed magnetic field. The 
gamma-ray luminosity due to pion decay is ~ 1.0 1 23 erg 
s _1 and the corresponding to relativistic bremsstrahlung, 
~ 4.7 1 29 erg s _1 . 

The EGRET threshold towards the position of WR 147 
is similar to that for WR 146. A gamma-ray source at the 
distance of the binary system would have been detected 
by EGRET if its luminosity is greater than 2.2 10 34 erg 
s -1 , i.e more than a factor two above the expected level of 
WR 147, according to the estimates based on the param- 
eters in Table |[ As for WR 146, this fact could explain 
why no gamma-ray source was found yet towards WR 147. 



Star 


B* 






CGI 


(G) 


WR 140 


440 


230 


WR 146 


1200 


30 


WR 147 


50 


4 



6. Discussion 

6.1. Stellar magnetic fields and the Razin-Tsytovich 
effect 



The magnetic field on the star surface can be computed 
using Eq. (|l]) and the magnetic fields estimated for the 
particle acceleration regions in Sect. I 

The rotational velocity is not well determined for the 
binaries selected here. We consider typical values of 400 
km s _1 for the WR stars, and r& ~ R*. The resulting 
values for the surface stellar fields are given in Table || 
as B*. They are between ~ 50 and 1200 G, in accordance 
with typical values estimated by Maheswaran & Cassinclli 
( |1994[ ) for WR stars. 

If the magnetic field is much stronger, then WR 140 
would fall below EGRET sensitivity. On the other hand, 
as indicated by White & Chen (1995), if the field is much 
weaker, the Razin-Tsytovich effect would suppress the ob- 
served synchrotron emission at cm wavelengths. The fre- 
quency at which the radiation is suppressed is: 



20 Hz, 



(12) 



where n c is the electron density expressed in cm~ 3 . 
This density is determined, in turn, by the stellar winds 
through: 



Mi 



47tr 2 v 



oo,lM m p 



(13) 



if rn p is the proton mass. 

We have taken the following minimum radio frequen- 
cies at which the binary systems have been observed: 1.4 
GHz for WR 140 (Wh ite fc Becker |l995| ), 327 MHz for 
WR 146 (Taylor et a l. ^996) , and 350 MHz for WR 147 
(Setia Gunawan et al. |2001a ). The derived minimum mag- 
netic fields near the shock front using Eq. ( |T^ ) are 100 mG 
for WR 140, 0.6 mG for WR 146, and 0.4 mG for WR 
147. The corresponding surface magnetic fields are given 
in Table || as -B™ 111 . These lower limits are quite consistent 
with our estimates. 

During the EGRET observations towards the WR 140 
system, it occurred a periastron passage, on 1993.2. 
The corresponding viewing period (212.0) lasted from 
March 9 to March 23, 1993, and the measured flux was 
25.7±12.7 10~ 8 ph cm" 2 s -1 . If we assume that all this 
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Table 5. Gamma-ray production in the colliding wind regions of WR 140, WR 146, and WR 147, for the energy range 
100 MeV < E < 20 GeV 



Stellar system 


Mechanism 


Expected 
luminosity 
(erg s" 1 ) 


Observed luminosity 
or upper limit 1 
(erg s" 1 ) 


WR 140 


Inverse Compton scattering 
Relativistic bremsstrahlurig 
Neutral pion decay 


~ 2.1 10 34 
~ 1.3 10 31 
~ 1.2 10 21 


~ 3.2 10 34 


WR 146 


Inverse Compton scattering 
Relativistic bremsstrahlung 
Neutral pion decay 


- 2.8 10 33 
~ 4.6 10 29 
~ 2.5 10 21 


< 8.0 10 34 


WR 147 


Inverse Compton scattering 
Relativistic bremsstrahlung 
Neutral pion decay 


~ 8.0 10 33 
~ 4.7 10 29 
~ 1.0 10 23 


< 2.2 10 34 



1 From EGRET observations and assuming the gamma-ray sources at the same 
distances of the stars. 



flux is due to IC scattering at the colliding wind region, 
the value of the magnetic field at the shock region must 
be about 0.16 G, which leads to a stellar surface magnetic 
field of ~ 350 G. 



6.2. Variability 

The size of the non-thermal regions in the WR binaries 



15 



cm, 



considered in this research is in the range 10 — 10 
which is about 3 orders of magnitude larger than the size 
of the stars. Hence occultation events cannot produce sig- 
nificant flux changes. Variability, instead, can result be- 
cause of the changing UV photon flux originated in very 
eccentric orbits and also from changes in the wind outflow. 
The timescales associated with the first process, however, 
exceeds the total EGRET observing lifetime by a factor 
of 2, at least in WR 140, till now the only detectable case. 
Future GLAST studies might reveal a modulation of the 
gamma-ray emission with the orbital period for this sys- 
tem. Changes in the winds could affect the injection and 
acceleration rates at the shocks on shorter timescales, but 
these effects are below the sensitivity of EGRET in the 
case of WR 140. Long term monitoring with high- altitude 
grou nd-bas ed GeV Cherenkov arrays like 5@5 (Aharonian 
et al. 2001 ) could be very useful to establish the changing 
properties of the winds. 

6.3. Additional gamma-ray contributions in the region 
of WR 140 

The interstellar matter surrounding the WR 140 system 
has been studied by means of HI-21cm radio and IR obser- 



vations (Arnal 2001 and references therein). Arnal found 
a minimum in the neutral hydrogen distribution, built by 



the action of the stellar winds, and estimated an HI mass 
of the surrounding emission in ~ 1300 solar masses. The 
HI void has a major axis of ~ 11.5 pc and a minor axis 
of ~ 8.4 pc. If particle re-acceleration is occurring at the 
terminal shock of the wind, then the material accumu- 
lated in the shell of the HI bubble could be exposed to 
relativistic proton bombarding yielding an additional con- 
tribution to the the total gamma-ray flux measured from 
3EG J2022+4317. This is a particularly interesting possi- 
bility, since the source is classified as "possibly extended 
source" in the 3EG catalog. 

The 7-ray flux expected from p — p interactions is 
(Aharonian & Atoyan 1996): 



F(E > 100 MeV) = — q 7 M m d' 

47T 



(14) 



where <7 7 is the gamma-ray emissivity of the medium 
whose total mass is Mm, d is the distance to the source, 



and 



'li- 



as before, is the proton mass. The gamma-ray 



emissivity can be related to the value observed in the 
vicinity of Earth by q y = kqQ. If the proton spectrum 
near the Earth is similar to the spectrum in the sur- 
roundings of WR 140 we can approximate k ss fccR, i.e. 
the enhancement factor is given by the ratio of cosmic 
ray densities in the shell and in the solar neighborhood. 
The latter density is usually taken as ~ 1 eV cm' 3 , and 
then q Q (E > 100 MeV) « 2.2 x 10' 25 (H- atoms)' 1 
(e.g. Dermer 1986| ). If we consider as reasonable a cosmic 
ray enhancement in the range k ~ 10 — 100 due to re- 
acceleration at the wind terminal shock and we adopt the 
estimate Mhi ~ 1300 Mq for the available mass (Arnal 
pooT| ), we obtain an additional contribution to the total 
gamma-ray flux from 3EG J2022+4317 of - 1.7-17 10' 8 
ph cm' 2 s . This extended emission could be responsible 
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for the confusion in the EGRET detection. This flux will 
contribute with a luminosity of ~ 5.7 10 32 - 10 33 erg s~ 4 . 

6.4. Remarks on the electron spectrum of WR 147 



Table 7. Gamma-ray flux expected at other energy ranges 



As we have mentioned in Section 3.3 it is not clear that the 
electron injection spectrum for this system can be repre- 
sented by a canonical value p — 2 as we have assumed for 
our calculations. Even taken into account the effect of the 
Compton losses it is not possible to explain the significant 
steepening observed in the synchrotron flux distribution 
above 10 GHz (Skinner et al. 1999 ). The radio spectrum 
is similar to what would be expected from a monoenergetic 
electron population. In this case, the spectra l distr ibution 
of the synchrotron radiation is (e.g. Longair 1997): 



#5/3 



(15) 



where K 5 / 3 (r]) is the modified Bessel function and v c is 
the critical frequency near which the emission has its max- 
imum. In the case that the monoenergetic electrons inter- 
act with an external photon field, the spectral distribution 
of the emerging IC radiation will mimic the seed spec- 
trum. With the photon fields considered in this paper, the 
expected gamma-ray spectrum will be similar to a black 
body spectrum. In this way, gamma-ray observations of 
WR 147 with good sensitivity and high spectral resolution 
can be used to probe the nature of the injected relativis- 
tic particle population and to test Skinner et al.'s ( 1999 ) 
proposal of a monoenergetic electron spectrum. Hopefully, 
GLAST will be able solve this problem. 

6.5. Predictions for upcoming satellites 

The gamma-ray fluxes expected at other energy ranges 
than EGRET's, due to IC scattering of UV photons from 
the secondary, are given in Table 6. It can be seen that 
some INTEGRAL and GLAST detections can be ex- 
pected. The continuum sensitivity of INTEGRAL'S IBIS 
instrument is 2 1CP 7 ph cm -2 s^ 1 at 1 MeV for an expo- 
sure of 10 6 s. GLAST sensitivity at E > 100 MeV for one 
year survey is ~ 4 10 -9 ph cm -2 s _1 . High-quality data 
from these instruments will help to fix the spectral shape 
of the sources in different energy ranges. In particular, no- 
tice that IBIS might observe the spectral break predicted 
at - 2.4 MeV for WR 140, whereas GLAST should detect 
the steepening in the spectrum of WR 147 at energies > 
70 GeV. 

7. Conclusions 

In this paper we have considered the production of 
gamma-ray emission in three WR+OB binaries that are 
well-known non-thermal radio sources. These systems dis- 
play clear evidence of a colliding wind zone where strong 
shocks are formed. The existence of a significant syn- 
chrotron radio emission from these regions ensures the 



System 


F 7 [15 keV -10 MeV] 1 
(ph cm -2 s -1 ) 


F 7 [20 MeV - 300 GeV] 2 
(ph cm~ 2 s -1 ) 


WR 140 


8.3 10" 4 


l.lxlO -6 


WR 146 


1.2 10 -4 


8.9xl0~ 8 


WR 147 


1.2 10~ 3 


9.4xl0~ 7 



1: INTEGRAL energy range; 2: GLAST energy range 



presence of locally accelerated relativistic electrons, and 
since the regions are also exposed to strong stellar pho- 
ton fields, the necessary conditions for inverse Compton 
production of high-energy photons are fulfilled. We have 
studied here whereas, according to the available multi- 
frequency information, this high-energy emission is de- 
tectable with current technology. We conclude that in the 
case of WR 140, the expected flux is strong enough as 
to account for an already observed but yet unidentified 
EGRET source: 3EG J2022+4317. In the case of the WR 
stars WR 146 and WR 147 the fluxes are below the current 
detection thresholds, but forthcoming experiments could 
detect them. 

If the existence of gamma-ray emission from early-type 
stars can be established beyond all observational doubt, 
the implications for galactic cosmic ray astrophysics could 
be very important. An important pending issue is what 
are the maximum energies at which protons could be ac- 
celerated by these systems. If re- acceleration takes place 
efficiently at the terminal shock fronts of the winds, then 
energies above 100 TeV could be expected. The recent de- 
tection of Cyg OB2 association by HEGRA instrument 
(Aharonian et al. 2002) could be a first step towards the 
identification of some stellar systems as cosmic ray sources 
below the so-called knee in the galactic cosmic ray spec- 
trum. 
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